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Using the published KTeV samples of Kl — » -K^e^u and Kl — > Tv^fi^v decays [l|], we perform 
a reanalysis of the scalar and vector form factors based on the dispersive parameterization 0, 01- 
We obtain phase space integrals 7|: = 0.15446 ± 0.00025 and I^ = 0.10219 ± 0.00025. For the 
scalar form factor parameterization, the only free parameter is the normalized form factor value at 
the Callan-Treiman point (C); our best fit results in InC = 0.1915 ± 0.0122. We also study the 
sensitivity of C to different parametrizations of the vector form factor. The results for the phase 
space integrals and C are then used to make tests of the Standard Model. Finally, we compare our 
results with lattice QCD calculations of Fk/Ftt and /+(0). 
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I. INTRODUCTION 



Recently, much effort has been devoted to measure the 
vector and scalar Kit form factors in semileptonic kaon 
decays in order to determine the phase space integrals. 
These integrals, along with the kaon branching fractions, 
allow to determine the GKM matrix element \Vus\- The 
scalar form factor (/o) is difficult to measure because it is 
kinematically suppressed in K^s decays, and is therefore 



only measurable in K^^ decays, which have contributions 
from both the scalar and vector form factors. In addi- 
tion, with the present experimental precision, only one 
fn model parameter can be accurately measured (see e.g. 
[^). Until recently, these form factors were determined 
with Taylor and pole parametrizations [l|, [^, d, 0| : 
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where t =■ [px ^ PttY is the expansion parameter as 
a function of the kaon and pion four-momenta, and 
f+,o{t) = /+.o(0//+-o(0) are the normalized vector and 
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scalar form factors. 



to the data are q and A" 



The parameters measured in a fit 
the slope of the form fac- 
tors and the curvature of the vector one, and My.s- the 
mass of the vector and scalar resonances. While the sec- 
ond order Taylor expansion has been used to measure 
the vector form factor with sufficient precision, the scalar 
form factor can only be determined using the first-order 
Taylor expansion or the pole model. However, clearly 
one has at least to know the curvature to have a proper 
description of fo{t) in the physical region of if^a-decays. 
After results were reported based on the Taylor and pole 
parametrizations, a form factor parametrization based 
on conformal mapping was discussed in Ref. Q in the 
context of i? ^ irli/ to improve the convergence of the se- 
ries and to give rigorous bounds on its coefficients. This 
parametrization was applied to the K^^ case in Ref. 
and recently used by the KTeV collaboration to reanalyze 
their K^s data [Toj . 

As an alternative approach, the dispersive 
parametrization in Refs. 0, Q has the advantage 
to account for the correlation between the slope and 
the curvature, by using low energy Kit phase shifts 
[HI, [13, Ell • It involves only one free parameter for both 
the scalar and vector form factors to be determined 
from the existing data sample. The sole scalar form 
factor parameter is C, the value of the normalized 
scalar form factor at the Callan-Treiman (CT) point, 
t = Aktt = rnj^ — iri^, the difference of kaon and pion 
masses squared. Once C is determined, the shape of the 
scalar form factor is known with a high precision in the 
physical region and somewhat beyond. The choice of 
this particular parameter C is guided by the existence 
of the Callan-Treiman theorem [iJl which predicts its 
value in the SU(2) x SU(2) chiral hmit. For physical 
quark masses. 



C = /o(A/i-7r) — 
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where F^^ and Fx are the pion and kaon decay con- 
stants, respectively, and Act is a correction of or- 
der O {mu 4/4,1: Ftt) arising from non zero quark masses 
nimmd- This correction has been evaluated within Chi- 
ral Perturbation Theory (ChPT) and is small enough 
that the right-hand side of Eq. ^ can be determined 
with sufficient accuracy as discussed in § IIVI to compare 
with C measured in iiT^a-decays. Thus apart from the de- 
termination of \Vus\, which is used to test the unitarity 
of the CKM matrix within the Standard Model (SM), 
a measurement of the scalar form factor at the Callan 
Treiman point provides another interesting test of the 
SM, namely a test of the couplings of light quarks to 
W . Another interest in the experimental determination 
of the shape of the Ktt scalar form factor is the possibil- 
ity of determining low energy constants which appear in 
ChPT [Hi. 

The NA48 [H and KLOE [Hi collaborations have re- 
analyzed their data with the dispersive parameterization 
[3l. The values of C obtained in these two experiments 



differ by 2.1a. Here we present a similar reanalysis of 
the KTeV data [l| leading to an improvement on the 
precision on the determination of the form factors com- 
pared with the previous KTeV results [l|, E3l . Since the 
vector and scalar form factors are correlated, alternative 
parametrizations for the vector form factor are studied 
to probe the robustness of the scalar form factor result. 

The paper is organized as follows. In § |TT] we present 
the results of the dispersive analysis of the KTeV data. In 
§ mil we discuss the correlations between the vector and 
the scalar form factor. ? IIVI is devoted to a discussion of 
different applications of our results, in particular the test 
of the SM. We summarize in S IVl 



II. DISPERSIVE ANALYSIS OF KTEV 
SEMILEPTONIC DATA 

Assuming fo{t) is never equal to zero, the dispersive 
representation for the normalized scalar form factor reads 



fo{t) = exp 
G(t) = 
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Ak.. 
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Note that C is here the only free parameter. (/>o(s) repre- 
sents the phase of the form factor: following Watson's 
theorem [14], this phase is equal to the Ktt scatter- 
ing phase within the elastic region. In writing Eq. ((4|. 
two subtractions have been made to minimize the un- 
known high energy contribution to the dispersive inte- 
gral, Eq. ([5]). The two subtraction points have been taken 
at t = and at the CT point to take advantage of the CT 
theorem, Eq. ([3]). The resulting fimction G{t) in Eq. ([5]) 
does not exceed 20% of the expected value of InC; since 
theoretical uncertainties on G{t) are ^ 10% its value, the 
corresponding uncertainty on InC is then a few percent 
of its value. 

The dispersive representation of the vector form fac- 
tor is constructed in a similar manner. Since there is no 
analog of the CT theorem in this case, the two subtrac- 
tions are performed at i = 0. The normalized vector form 
factor is 
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where A-|_ = irL^df^[t)/dt\t=o and (j)+{s) is the phase 
of the vector form factor. As in the case for the scalar 
form factor, information on the Ktt phase shifts in the 
elastic region is used to determine 0_|_(s). The main con- 
tribution to (t>+{s) is the dominant /C*(892) resonance. 
The extrapolation of the Kn phase shift data down to 
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Kez only 


K^-i only 


Ae3 and K^'i Combined 


A+ X 10^ 
InC 

p(A+,lnC) 
/dof 


25.17 ± 0.62 
66.6/65 


24.57 ± 1.10 
0.1947 ± 0.0140 
-0.557 
193/236 


25.09 ± 0.55 
0.1915 ± 0.0122 
-0.269 
0.48/2 


AV X 10' 

A_|_ X lU 

Ao X 10' 


25.17 ± 0.62 

i.ZZ zt U.U ( 


24.57 ± 1.10 

l.ly =h U.U / 

1 q 99 _|_ 1 on 

0.59 ± 0.05 


25.09 ± 0.55 
i.zi zt u.Uo 

1 9 _|_ 1 17 
iz.yo Jz 1.1/ 

0.58 ± 0.05 




0.15450 ± 0.00028 


0.15416 ± 0.00060 
0.10207 ± 0.00032 


0.15446 ± 0.00025 
0.10219 ± 0.00025 


I'k/Ik 




0.6621 ± 0.0018 


0.6616 ± 0.0015 



TABLE I: Results of the analysis of the KTeV KlsS and Kl 113 data using a dispersive parameterization for the vector and scalar 
form factors. A+ and In C are the parameters of the fit used to calculate the slopes, curvatures and the phase space integrals. 
The uncertainties correspond to the total ones, adding the statistic, the systematic as well as the theoretical ones in quadrature, 
see the text for more details. 



threshold is done here following a Gounaris-Sakurai con- 
struction based on the K*{892) and exhibiting the cor- 
rect threshold behavior and the correct properties of an- 
alyticity and unitarity. The value of H(t) represents at 
most 20% of the value of A_|_ such that the latter can be 
measured with high precision. For more details on the 
dispersive representations, see Refs. and 0]. 

In Rcf. 3, a thorough discussion of the different 
sources of theoretical uncertainties of the dispersive rep- 
resentations can be found. They include the error on the 
low energy Kt: phase shifts and an estimate of the uncer- 
tainties due to the unknown high energy behaviour of the 
phases 0o(s) and (s). The corresponding error-bands. 
6G{t) and 6H{t), are used in this analysis to propagate 
uncertainties on A_|_ and In C. 

The analysis of the KTeV data is done using their KlcS 
and i^'i^s samples with 1.9 x 10^ and 1.5 x 10^ events, 
respectively after selection requirements. These samples 
were collected in a special run in which the beam in- 
tensity was lowered by a factor of ten compared to that 
used to measure e'/e. The laboratory-frame kaon ener- 
gies are 40-160 GeV (mean is 70 GeV), and the momenta 
of charged particles are measured with much better than 
1% precision. Muons are identified with a large scintil- 
lator hodoscope behind 3 meters of steel. Electrons and 
pions are identified primarily by ratio of energy deposited 
in the cesium iodide calorimeter (E) to the momentum 
measured in a magnetic spectrometer (p); E/p 1 for 
electrons, and E/p < 1 for pions. In addition to us- 
ing the KTeV data, we also use the KTeV Monte Carlo 
(MG) to correct for the detector acceptance that results 
in a non-uniform sampling of the Ki^ Dalitz plot. 

The results of the dispersive analysis are given in Table 
m The associated slope and curvature are also given, 
based on Taylor expansions of Eqs. ^ and © using the 
best- fit values of In C and A+ , respectively. A combined 



Ke3 and Jf^3 dispersive analysis leads to 

A+ = 0.02509 ± 0.00035stat ± 0.00027syst ± 0.00033th 
= 0.02509 ±0.00055, 

InC = 0.1915 ±0.0078stat±0.0086syst± 0.0038th 
= 0.1915 ±0.0122. 

(8) 

The phase space integrals are then given by: 

11^ = 0.15446 ±0.00019stat± 0.00015syst± 0.00008th 

= 0.15446 ±0.00025, 
I^^ = 0.10219 ±0.00017stat±0.00017syst± 0.00005th 

= 0.10219 ±0.00025. 

. (9) 

Systematic uncertainties are estimated following the pre- 
scription of Ref. [l3| , by scaling the ratio of systcmatic- 
to-statistical uncertainties for the pole model in Table 1 
of Rcf. Q (seeEq. (15) of Ref. Statistical, system- 

atic and theoretical uncertainties are added in quadra- 
ture to give the total uncertainty. To estimate the the- 
oretical error on In C and A+ induced by uncertainties 
on the functions G{t) and H{t) entering the dispersive 
representations, we perform fits using Git) ± SG(t) and 
H{t) ± 6H{t). The function G{t) is positively correlated 
with InC, and H{t) is negatively correlated with A_|_; 
these correlations lead to reduced uncertainties in the 
phase-space integrals {Ik)- 

Table U provides also values of the phase space inte- 
grals ratio, /^//J^. Note that for them the estimated 
total uncertainty takes into account correlation due to 
the common vector form factor f+{t) which reduces the 
uncertainty 

After subtracting the common theoretical uncertain- 
ties, our result for InC is consistent with the KLOE 
resuh, InC = 0.2038(246) [l3|, and it is 2.6cr larger 
than the NA48 result, InC = 0.1438(140) [ll|. For 
the previous form factor fits from KTeV [l|], the phase- 
space integrals (/^ = 0.15350(105),/^ = 0.10165(80)) 
arc in good agreement with the dispersive results above, 
but the precision is limited by modeling uncertainties 
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that are twice as large as the statistical errors. The 
large modeling uncertainty in Ref. [l[ is based on 
the difference between using the pole model (/J^ = 
0.15U5{23stat),lK = 0.10235(22stat)) and the quadratic 
model (/^ = 0.15350{Ustat),lK = 0.10165(39,tat)), 
where the uncertainties are statistical only; the pole 
model result agrees very well with the result above based 
on the dispersive analysis. For the z-parameterization 
[3 [Ik 0.15392(48)), the theoretical uncertainty is 
slightly larger than that from the dispersive analysis, 
and therefore the /^-discrepancy may be significant. To 
investigate this difference, several Monte Carlo samples 
were generated using input form factors from the result 
of the dispersive fit and subsequently analysed using z- 
parameterization. Based on this study, the Ik integrals 
obtained with the z-parameterization reproduce on aver- 
age the input value, and the difference between the z and 
dispersive parameterization observed for the KTeV data 
is consistent with a l.SCTgtat fluctuation. 



III. COMPARISON OF DIFFERENT FORM 
FACTOR PARAMETERIZATIONS 

As pointed out in the introduction, a main advantage 
of the dispersive parameterization is the possibility of 
determining the value of the scalar form factor at the 
Callan-Treiman point, thus allowing for a test of the SM. 
At present, only the dispersive parameterization makes it 
possible to determine the scalar form factor at this point, 
which lies far beyond the cndpoint of the physical region, 
with reasonable precision. Since there is a large correla- 
tion between In C and the slope of the vector form factor 
(see p( A_|_ , In C) in Table |T| , however, it is important to 
study the sensitivity of In C to the choice of parameteri- 
zation for the vector form factor. 

To investigate this sensitivity, we have fit the KTeV 
data using the dispersive parameterization for the scalar 
form factor (with In C as a free parameter) and four dif- 
ferent parameterizations for the vector form factor: 

• the dispersive parameterization Eq. ^ 

• the pole parameterization Eq. ^ 

• the quadratic (second-order) Taylor expansion 
Eq. KID 



the z-parameterization 



F+it)^F+ito) 



(^+(to,to,Q^) 

4>+ it,to, 



fe=0 

f+{t)^F+it)/F+iO) , 



(10) 



based on a conformal mapping of t onto the variable 
z with 



(11) 



In Eq. ([TOl) we have used the notation from Ref 
i. 

The results of these fits are summarized in Tables [IT] 
and mil The "dispersive" results are taken from Sec- 
tion ini All of the fits have good /dof. Interest- 
ingly, the pole and dispersive parameterizations result 
in very similar values for InC, while the quadratic and 
z-parameterization results are similar. This can in fact 
be easily understood from the Appendix which presents a 
detailed investigation of the correlations between param- 
eters in the different parameterizations. All of the InC 
results are consistent within 2(Tstat, as can be deduced 
from the difference of the (4 units per one degree of 
freedom change) and by estimating the uncertainty of the 
difference as a difference of the uncertainties in quadra- 
ture: In C|disp/disp - In C|2/disp = 0.026 ± 0.013. 

A similar level of agreement between the different 
parametrizations is observed for the Ik integrals. For 
example for K^^^ data, using the z-parameterization 
with two parameters for each of the vector and scalar 
form factors, one obtains Ik\z = 0.15331 ± 0.00072stat ± 
0.00040syst and I^\,^, = 0.10101 ± 0.00053stat ± 
0.00039syst. 



IV. DISCUSSION 

While both the z and dispersive parametrizations give 
rigorous boimds on theoretical uncertainties, the latter 
uses additional experimental input such as the low en- 
ergy Kn phase shifts. This allows for a one-parameter fit 
of the vector and scalar form factors, resulting in smaller 
uncertainties. In the following, we discuss the impact of 
the value obtained for In C, in the dispersive parametriza- 
tion. 



A. Comparison with lattice results 

Here we compare our result for In C against the lattice 
QCD calculations. This comparison does not depend on 
SM assumptions since no electrowcak couplings are in- 
volved. 

Figure [1] shows lattice QCD results for Fk/F„ and 
/+(0) in the 2+1 flavor case [H, [H, HO, [SI, [^, [H . A 
first classification of these data can be found in the recent 
proceeding Ref. [l^l awaiting for the FLAVIAnet Lattice 
Averaging Group's one. We have only considered pub- 
lished results and showed them as bands including sys- 
tematic and statistical errors not giving the central values 
for clarity. Note that RBC/UKQCD has much bigger 
systematic errors compared to the other collaborations 
leading to the rather large band for Fk/F^^ ranging from 
1.14 to 1.27. 

Also shown is the /+(0) vs Fk/F.„: dependence as de- 
rived from Eq. ([3]) using our result for In Cj42| . For Act, 



5 





Parameterization Vector FF/ Scalar FF 


Results 


disp.(I)/disp. 


pole(I) /disp. 


quad (II)/disp. 


z-param.(II) /disp. 


i. Ill 1JCU.C1111, t/j 


A+ = 24.57(83) 


Mv = 890.00 


\+ = 17.5(3.4) 


ai = 1.057(63) 






(13.00)MeV 


A+ =4.3(1.4) 


aa = 3.9(3.2) 


A', 


24.57(83) 


24.59(72) 


17.5(3.4) 


20.00(2.60) 


II 


1.19(4) 


1.21(7) 


4.3(1.4) 


2.5(6) 


Fit param. In C 


0.1947(91) 


0.1944(93) 


0.169(16) 


0.170(16) 


Ao 


13.22(78) 


13.20(79) 


11.03(1.37) 


11.11(1.37) 


a; 


0.59(2) 


0.59(2) 


0.54(3) 


0.54(3) 


p{vi, In C) 


-0.557 


0.588 


0.707 


0.477 








-0.819 


-0.766 


x7dof 


193/236 


193/236 


189/235 


189/235 



TABLE 11: Results of the fit to Kni3 data using different parameterizations for the vector form factor and the dispersive one 
for the scalar form factor. The uncertainties are only statistical. A+, A+, A", Ao and Aq are in units 10^. 





Parameterization Vector FF 


Results 


dispersive (I) 


pole (I) 


quadratic (II) 


z-param. (II) 


Fit param. 


A+ = 25.17(38) 


Mv = 881.03 


X'+ = 21.67(1.59) 


ai = 1.023(28) 






(5.91)MeV 


A+ = 2.87(66) 


aa = 0.75(1.58) 


AV 


25.17(38) 


25.10(41) 


21.67(1.59) 


22.69(1.20) 


A'l 


1.22(4) 


1.26(3) 


2.87(66) 


1.93(30) 


Correlation 






-0.96 


-0.064 


xVdof 


66.6/65 


66.3/65 


62.2/64 


62.3/64 



TABLE III: Results of the fit to Kles, data using different parameterizations with only the statistical uncertainties. Pole, 
quadratic and z -parameterization results are from Refs. Q/ and JTd^ . A+, A^ and A+ are in units 10''. 



we have used the value |2£ 



Treiman relation, now rewritten as 



Act = -0.0035 ±0.0080, 



(12) 



taken from a next-to-leading-order calculation in chiral 
perturbation theory in the isospin limit. The error is a 
conservative estimate of higher order corrections in the 
quark masses 77i„, rrid and TOs [2^. This value of Act 
is in agreement with other recent determinations [2^ . [28l . 

in,!!!. 

Combining all the lattice results, the grey band in 
the ./+(0) vs Fk/Ftj plane shown in Fig. [1] is obtained. 
Comparing this band to the KTeV result, the ranges 
Fk/F^ < 1.20 and /+(0) > 0.96 are favoured by the 
KTeV data. 



Fk 



= C- A 



CT- 



(13) 



This test consists of comparing the value of 
Fk/{Ftt ■ f+{0)), deduced from the Klis disper- 
sive form factor parametrization fit, to the value of 
Fk /{Ftt ■ f+{0))\sM , determined by assuming the SM 
electrowcak (EW) couplings and using the experimental 
(photon inclusive) branching fractions /^-n^ ^-^^ 
Tkl^3 measurements. 
We define r as 



(C - Act) 



F.-f+{0) 
Fk 



SM 



(14) 



B. Test of the SM 

As mentioned in § [H the small size of the Act cor- 
rection allows for an accurate SM test using the Callan 



Physics beyond the SM, such as modifications of EW 
couplings of quarks due to new exchanges close to the 
TeV scale, could cause r to differ from unity. An example 
of modified EW couplings between right-handed quarks 
and the W boson is discussed in Refs. 0, [3l|. 

We first calculate Fk/{Ftj ■ /+(0))|sm- Assuming the 
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1.02 - 



0.98 - I 



0.94 



0.92 



T 



T 



KTeV 

HPQCD-UKQCD 

NPLQCD 

MILC 

RBC-UKQCD 



I 



0.96 - v-'^e- - 



t; 



L 



1.14 



1.16 



_L1j_ 



1.18 



1.2 
F^/F 

K Jl 



1.22 



1.24 



1.26 



FIG. 1: The constraints on /+(0) and Fk/Ft^ from vari- 
ous lattice QCD calculations. The KTeV result for Fk:/(-FV • 
/+(0)) is shown as the red solid line, derived by using the 
Callan Treiman relation and the NLO result for Act- Also 
shown is the error due to Act (red dotted lines), and the re- 
sulting error when added in quadrature to the total error on 
C as given in Eq. 0) (red dot-dashed line). 



SM couplings, one has 



\SM 



r[7r Ivii'^)] 

\V.us? FlmK{l-xlY [1 + ^F{xk)] 



cte- 



(15) 



where xp = mi /Alp. The expression for cte, which de- 
pends on the hadronic structure and particle masses, can 
be found in Ref. [s^. The function F{x) parametrizes 
the electromagnetic radiative corrections, and a is the 
fine structure constant. The coefficient A4 thus defined 
is equal to 0.2387(4) (see Ref. 

The K^es partial width is expressed as 



Kl 



\sM = M^\f+{0)V,s\' 



where 



« /-2 _ ^2 

^ "^^(192^ 



Sew (1 + S'k) r 



K- 



(16) 



(17) 



Here Gp is the Fermi constant, Spw are short-distance 
electroweak corrections, and denotes the electromag- 
netic (EM) radiative corrections. From equations ([1| 
and (I16D. it can be shown that 



Fk 



F.-f+iO) 



r 



1/2 



SM 



r + -r 



\Vud\ 



Kl 



M' 



(18) 



Using the world average result [33| /r^4 
1.3337(46), the KTeV measurement of Tk^^s 



0.4067(11) Q, the values of from Table U and 5 



0.0130(30) from Ref. [39|, and the value of |Kd| inferred 
from 0+ 0+ superallowed nuclear transitions [13] [43 



iKdl = 0.97418(26), 
Eq. (flS)) gives the result 
Fk 



F.-f+iO) 



1.2407 ±0.0044. 



SM 



(19) 



(20) 



This result can be compared with the experimental 
determination of FxliF^^ ■ /+(0)) through K^a decay as 
given in Table 1. One obtains 



r = 1.0216 ± 0.0124exp ± 0.0039thco ± 0.0067ac 



(21) 



The first two errors come from the experimental and 
the theoretical uncertainties on InC respectively, and 
the last error comes from the estimated error on Act 
[44| . Adding the different errors in quadrature, we ob- 
tain r = 1.0216 ±0.0146. 

Analogous to Fig.[l] Fig. [2] shows the two bands in the 
f+{Q)-vs-FK / Ft^ plane. The first band (red) shows the 
dispersive parametrization analysis of the KTeV 
data. The second band (green) shows the SM predic- 
tion, Eq. (pn)l . We observe a 1.5 a difference between the 
KTeV result and the SM prediction. 



0.95 




1.15 



1.25 



FIG. 2: The constraints on /+(0) and Fk/Ft^. The inserted 
— Xaiin ~ 1 ellipse is the SM result obtained from Eqs. 
f_?5|) and \16\] and the CKM unitarity. The dashed green line 
and its error bar (green dotted line) corresponds to FkHFtt ■ 
/+(0))|sj\/, calculated using the measured decay widths and 
assuming the SM couplings of light quarks, Eq. f20fl . Shown 
in red is the Fk /{Ft^- f+{0)) result, derived from the dispersive 
parameterization fit to KTeV data and the NLO result for 
Act. Also shown is the error due to Act (red dotted lines) 
added in quadrature to the total error on C as given in Eq. 
(red dot-dashed line). 



Separate bounds on Fk/FttIsm 
derived from CKM unitarity [3ll . 



and /+(0)|sM can be 
|40| . Unitary implies 
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that IKtdP + IKisP + IKibP = 1, and measurements in- 
volving b ^ u transitions have shown that is negh- 
gibly small. Consequently, the SM mixing of light quarks 
is entirely specified by the value of \Vud\- 

Substituting |KisP = 1 — \Vud\^ into equations ([15]) 
and (fT6|) . and solving for Fk/Ft^ and /+(0), we obtain 
the contour in Fig. [21 One has 

Fk/F^\sm - 1.189 ±0.007 
f+msM = 0.959 ± 0.006. (22) 

With the current experimental precision, the data show 
a marginal agreement with the SM as concluded before. 



The measured value of In C, the scalar form factor at 
the Callan-Treiman point, leads to a dependence of /+(0) 
on Fk/Ftj within 1.5(7 of the Standard Model prediction. 
It favors an Fk /^V value on the lower side of the lattice 
results, and an /+(0) value on the higher side. New val- 
ues of the decay phase space integrals and If^ are 
obtained, where the latter is consistent with the result 
obtained by z parametrization p^ . These new values 
can be used to determine the ratio of the decay constants 
of the two semileptonic modes, Gp/Gp, which is in excel- 
lent agreement with the Standard Model prediction. A 
detailed analysis of a Taylor-expansion vector form fac- 
tor fit to the data is used to study how the scalar and 
vector form factor correlations affect the result for In G. 



Ratio Cp/Gl 



Taking the ratio of the KlcS to Kj^^^ partial widths, 
without assuming the equality of the Gp'^ decay con- 
stants, one obtains 



L:rp 

G% 



r(ifi 



1 



Using the value of 
tion estimates (1 



k/Ik 



from Table 1 



iK_ 

(23) 

the EM correc- 



5|f) = 1.0058 ±0.0010 from 



Ref. |39|,|45| and the direct measurement of T{Kl — > 
Tr^l^'^iy)/r(KL 7r±e=Fi/) = 0.6640 ± 0.0014 ± 0.0022 
from Ref. [l|, one obtains 



7^ 



(jrp 

G% 



0.9978 ±0.0049. 



(24) 



This result is in excellent agreement with the Standard 
Model expectation of unity, and it is very similar to the 
previous KTeV result, {G^/Cpf = 0.9969 ± 0.0048. Al- 
though the I'^ and phase-space uncertainties in Ta- 
ble [Ij are smaller than in the previous KTeV analysis [l|, 
the uncertainty on TZ^/e is almost identical in these two 
analysis. This is due to the fact that the uncertainties on 
the radiative corrections dominate the uncertainties on 
this ratio. In both analyses, the uncertainty on the ratio 
of phase-space integrals, Ix/^k^ significantly smaller 
than the quadrature-sum of the individual uncertainties 
because of correlations in the vector form factor (/+). In 
the KTeV analysis, the theoretical uncertainty related to 
the scalar form factor (/p) could not be evaluated be- 
cause of the parameterization used. In this analysis us- 
ing the dispersive parametrization, the /o uncertainty is 
more reliable, resulting in a more robust estimate of the 
uncertainty on 7?.^/e. 



V. SUMMARY 

A dispersive analysis of the semileptonic form factors 
for /•^Les and Kl^s has been performed based on the 
published KTeV data. 



APPENDIX A: PARAMETER CORRELATIONS 

In this Appendix, we will investigate correlations be- 
tween parameters in the different form factor parameter- 
izations. This study helps to understand several results 
discussed in the text, especially the difference between 
the parametrizations used for the vector form factor pre- 
sented in Sec. IIIIl and the robustness of the dispersive 
result presented in Sec. [TTj 

For this study, we will perform several fits of the A'^^a 
data. We will always use the dispersive parameterization 
for the scalar form factor. For the vector form factor, 
we will consider a cubic expansion, i.e., the first three 
terms of Eq. ([T]) will be taken into account. Indeed, in 
the physical region of A'^^a decay, a good representation 
of the dispersive parameterization may be obtained by 
Taylor expanding it with respect to t and keeping only 
the first three terms. We thus have four parameters (In C, 
A^, A" , and A"') which will enter the fitting procedure. 

Let us first fix the third order coefficient A"' to zero, 
which corresponds to a fit with a quadratic parametriza- 
tion for the vector form factor. We then perform a fit 
with two free parameters, namely InC and A^, while 
fixing the curvature A" to three different values: A" = 
0.001, 0.002 and 0.003 (A^ = 0.0 would correspond to a 
linear parameterization for the vector form factor, which 
gives a very poor fit to the data) . The results of these fits 
are shown by the black points on Fig. [3| The value of A'_,_ 
can be read on the y axis and In G is given by the first 
number in parentheses above each black point; the second 
number in parentheses gives the of the fit. The linear 
relationship between A^ and A" seen on the figure shows 
the strong correlation between these two quantities. For 
example, a fit to K^s data using a quadratic parameter- 
ization for the vector form factor leads to a correlation, 
/o(A^,A") = —0.96 [l|. A similar hnear dependence is 
also obtained between InC and A'^_ or A" . Furthermore, 
when moving toward larger A" , In C decreases while the 
decreases only slightly. The green square at the end of 
the solid black line represents the minimum x^; it corre- 
sponds to the result of the quadratic/dispersive fit given 
in Section IIIIl (Wc have stopped the line at that value 
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0.02 



1 ' 1 

/ ' 

/ 

/ 


1 ' 1 ' 


# # quad./disp 






A- ■ -A cubic (A,^'"=0.U()08) /disp. 




/ 

/ 


■ fit quad./disp. 




/ 


^ fit disp./disp. 




/ 

/ 

/ 

(0.1^,194.7) 


relation slope/ciirvaure in disp. pat. 






[ — \ (0.190,191.9) 
(*.180,189.5p^«:^ 




/ (0.171, 18aOT^^^181, 190.1) 
MIN ■~----.a(0'166,18??)^^^ 


1 


^^^^0.1 69,1 8£ 
I.I. 


■2) 



I I I I I I I I I 

0.001 0.002 0.003 0.004 0.005 



V 

FIG. 3: Dependence of the fitted values of In C and A+ on the 
curvature, A" . The value of In C and the of the correspond- 
ing fit are given in parenthesis above each point. To simulate 
the presence of higher order terms in the z-parameterization, 
for the red line a third order term has been included. The 
effect of this third order term is only to shift the curve down- 
wards. In the legend, "quad./disp." for example indicates that 
a "quadratic" parametrization is used for the vector form fac- 
tor and a "dispersive" one for the scalar. 



of the curvature for convenience.) 

Next, we consider the impact of a third-order term 
in the Taylor expansion. Indeed, a dispersive or z- 
paramcterization of the vector form factor gives a non- 
zero value for the third order coefficient A"' when Tay- 
lor expanded. We repeat the same procedure as above 
(i.e., floating InC and A^, and varying A"), but the 
third-order coefficient is now fixed to A"' = 0.0008, 
close to what is obtained from the results of the z- 
parameterization (Tables II and III). The nonzero A"' 
causes an offset in the linear dependence of A'^ with A" , 
but leads to the same slope and a similar variation of In C 
with the slope and the curvature (see the red dot-dot- 
dashed line on Fig. [3]) [i^ . The minimum fit (green 
square for A"' = 0) is now moved toward a smaller A" and 
a larger A^ (red triangle denoted by 'MIN'. This point 
corresponds to the z-parameterization result in Tabic II.) 



The parameter shifts are required in order to compensate 
for the additional cubic term. Note that the value of in C 
is roughly the same at this minimum as the one obtained 
previously for A"' = 0. 

In the light of this study, let us now consider the other 
parameterizations used in the literature [l|, p^ . 
nj. One can distinguish two classes. One class (Class 

I) , of which the dispersive and pole parameterizations are 
examples, impose physically motivated relations between 
the slope and the curvature (and possibly the higher or- 
der terms in the Taylor expansion, the third order term 
being the most relevant one in the physical region). In- 
deed, as already emphasized, in the dispersive param- 
eterization the curvature and all higher order terms of 
the vector form factor are constrained not only on first 
principles, such as analyticity and unitarity, but also by 
including the information on the low energy Kn P-waves 
(K* resonance). The relation between slope and curva- 
ture is illustrated by the blue dashed curve in Fig. O 
This curve crosses the black solid line and the red dot- 
dot-dashed line for large values of A^ , small values of A" , 
and consequently large values for In C, nicely illustrating 
the result of the dispersive/dispersive fit (magenta dia- 
mond in the figure). In contrast, the second class (Class 

II) , of which the Taylor series and z-parameterization are 
examples, is based on mathematically rigorous expan- 
sions, in which the slope and curvature are free param- 
eters. Clearly, the existence of a relation between slope 
and curvature strongly constrains the fit in Class I, while 
the fit has more freedom when using Class II parameter- 
izations. One thus expects smaller fitting uncertainties 
in the results from the Class I and larger ones from Class 
II, so that if the theoretical errors are well controlled, the 
overall uncertainty will be smaller for class I. 
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Note that we have plotted /+(0) up to 1.02, though the 
Fubini Furlan bound claims /+(0) < 1. However, in 
QCD this bound can only be proven in the large Nc limit, 
leaving the possibility for a slight deviation from it |3^ . 
Note that a recent update [s^ has appeared. 
Note that this analysis of the KTeV data holds if there 
are no scalar, pseudoscalar and tensor couplings. Indeed 
in the presence of the two former what has been measured 
in table 1 is not In C but In C plus a small quantity, see 
Ref. [l^]. Consequently in that case a determination of 
r is not possible without a precise knowledge of these 
couplings. 

To maintain consistency of this analysis with that of m , 
the EM correction ratio has been taken from Ref. [391 ]. 
Recently, new radiative correction estimates have ap- 
peared based on ChPT calculations [JJl, which are in 
agreement with the old results but with a larger uncer- 
tainty. 

Note that for illustration purpose we only show a portion 
of this line as well as of the black solid line corresponding 
to the fits we have made. They of course go on in both 
directions. 



